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This  paper  investigates  two  proposed  methods  for  estimating 
the  condition  number  of  a matrix  from  factorizations  commonly  used 
to  solve  linear  systems.  One  method  estimates  the  condition  number 
with  respect  to  the  1-norm  from  the  LU  factorization,  and  the 
other  the  condition  number  with  respect  to  the  2-norm  from  the  QR 
factorization.  Random  matrices  of  various  orders  having  known 
distributions  of  singular  values  were  generated  and  the  estimated 
condition  numbers  compared  with  the  true  ones.  For  the  classes  of 
matrices  tested  in  this  study,  the  estimators  performed  rather  well, 
never  underestimating  the  condition  number  by  a factor  of  more  than 
ten.  The  paper  also  gives  an  efficient  method  for  generating  random 
^^arthogonal  matrices  with  the  Haar  distribution.  ^ ^ 
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Estimating  Condition 
Numbers  — an  Empirical 
Study 

G.  W.  Stewart 


1.  Introduction. 

This  paper  is  concerned  with  estimating  the  condition  number  of  a 
nonsingular  matrix  A of  order  n.  Let  ll-ll  denote  a vector  norm  and  its 
subordinate  matrix  norm[7_i.  The  condition  number  of  A is  the  number 

(1.1)  <p(A)  = ||  A llpll  A~  1ll  p . 

The  importance  of  the  condition  number  lies  in  the  following  theorem  (see 
[7,8]  for  proofs  of  a slight  variant). 

Theorem  1.1.  Let  x be  the  solution  of  the  system 

(1.2)  Ax  = b. 


Let  E be  a matrix  of  order  n satisfying 


‘p(e>  Pi?  < »• 

Then  A+E  is  nonsingular,  and  the  solution  of  the  system 

(A  + E)  x = b 

satisfies 


(1.3) 
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Theorem  1.1  allows  us  to  assess  the  accuracy  of  a solution  of  (1.2) 
when  the  matrix  A is  perturbed  by  the  addition  of  E.  Informally,  the 
bound  (1.3)  says  that  the  relative  error  in  A due  to  E is  magnified  by  a 
factor  of  <p(A)  in  the  solution. 

In  principal,  one  can  estimate  <p(A)  by  computing  A"1  and  esti- 
mating the  norms  in  (1.1).  This  approach  has  two  drawbacks.  First,  the  norms 
themselves  may  be  difficult  to  calculate.  More  important,  in  many  applications 
it  is  expensive  and  unnecessary  to  compute  A~*.  For  example,  the  usual  pro- 
cedure for  solving  linear  systems  such  as  (1.2)  is  to  factor  the  matrix  A 
into  the  product  of  two  computationally  more  tractible  matrices  and  use  this 
factorization  to  solve  the  system. 

One  way  of  circumventing  the  difficulties  raised  in  the  last  para- 
graph is  to  attempt  to  estimate  the  condition  number  of  A from  one  of  its 
factorizations.  The  purpose  of  this  paper  is  to  present  the  results  of  simu- 
lation studies  for  two  such  methods.  The  first  estimates  a condition  number 
from  the  LU  factorization,  the  second  from  the  QR  factorization.  The 
methods  are  applied  to  a variety  of  randomly  generated  matrices  with  known 
condition  numbers,  and  the  estimates  are  compared  with  the  true  values. 

Speaking  broadly,  we  conclude  that  the  methods  are  rather  good  for  the  class 
of  matrices  considered,  never  underestimating  the  condition  number  by  a 
factor  greater  than  ten. 

The  paper  is  organized  as  follows.  In  Sections  2 and  3 we  briefly 
describe  the  methods  under  consideration.  In  Section  4 we  describe  how  the 
random  matrices  used  in  the  study  are  generated;  in  particular,  we  give  an 


I 
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algorithm  for  generating  a random  orthogonal  matrix  that  is  in  some  sense 
uniformly  distributed  on  the  space  of  orthogonal  matrices.  In  Section  5 
we  describe  our  experiments  and  present  the  results  and  our  conclusions. 

2.  Estimation  of  k^. 

In  this  section  we  will  be  concerned  with  estimating  the  condition 
number  for  the  1-norm  defined  by 

llxilj  = i"=1  | xil. 

The  subordinate  matrix  norm  is  given  by 

/,  HAIL  = max  II  Axil,  = max  Z | a - . I . 

'Z#1'  1 ilxllj-1  1 l<j<n  i=l  J 

From  the  second  characterization  in  (2.1),  it  is  seen  that  ttA^  is  easily 

p 

computed  in  0(n  ) arithmetic  operations.  Thus  the  problem  of  estimating 
Kj  reduces  to  that  of  estimating  |lA~*llj. 

The  method  used  here  is  an  elaboration  of  one  described  in  l4J.  It 
has  been  treated  in  detail  elsewhere  [ 1 J , and  here  we  give  only  a brief  sketch. 
The  underlying  idea  is  that  the  maximum  in  the  middle  term  of  (2.1)  is  likely 
to  be  nearly  attained  for  any  vector  x.  Thus  to  approx  ,.jte 

||  A-1l|,  = max  II  A-1ull, , 

1 llullj-1  1 

we  start  with  a suitably  chosen  vector  u and  solve  the  system 


I 

i 


(2.2)  Av  = u, 

after  which  llA”1!!,  is  estimated  by  ||v||,/  Hull,. 
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Two  modifications  are  incorporated  into  this  basic  idea.  The  first 
insures  a high-quality  starting  vector  u.  If  Gaussian  elimination  is  used 
to  solve  the  system  (2.2),  the  matrix  A,  or  rather  a row  permutation  of  it, 
is  decomposed  in  the  form 


A = LI), 

where  L is  unit  lower  triangular  and  U is  upper  triangular.  The  system 
is  then  solved  by  succesively  solving 


Lu'  = u 


and 


Uv  = u'. 


Because  L is  lower  triangular,  the  element  u.j  depends  only  on  Uj,...,uj_j 
and  u.j.  Consequently,  the  choice  of  ui  can  be  deferred  until  ujj  has 

I 

been  computed,  at  which  point  it  may  be  chosen  to  enhance  the  size  of  u... 

The  second  modification  is  to  pass  once  more  through  the  matrix, 

solving 


.T 

A w = v, 

after  which  IIA'1!^  is  estimated  by  :| wli^/ilvil^,  where  iiwu^  = max  [ \ . 

The  reasons  for  this  step  are  detailed  in  (1);  essentially,  the  first  pass 
produces  a good  starting  vector  for  the  second  pass,  where  the  actual  estima- 
tion is  done. 
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A slight  variant  of  the  method  described  above  has  been  incorporated 
into  the  UNPACK  programs  for  solving  linear  systems  [2],  and  it  is  this 
version  that  we  study  in  this  paper. 

3.  Estimation  of 

In  principal  the  method  outlined  in  Section  2 can  be  used  to  estimate 

the  condition  number  with  respect  to  any  easily  evaluated  vector  norm  (an 

estimate  of  |lAllp  can  be  obtained  in  much  the  same  way  as  that  of  |IA~*!lp, 

if  a computationally  convenient  characterization  does  not  exist).  However, 

2 2 

for  the  2-norm  defined  by  ||xil2  = ex*,  there  is  an  alternative  method  based 
on  the  QR  factorization  with  column  pivoting.  Since  this  factorization  is 
often  computed  in  the  course  of  solving  least  squares  problems  (e.g.  see  L3, 
6,7  1),  we  include  a study  of  this  alternative  way  of  estimating 

A computational  method  based  on  Householder  transformations  for 
computing  the  QR  factorization  is  closely  connected  with  the  method  for 
generating  random  orthogonal  matrices  to  be  described  in  the  next  section. 
Consequently  we  sketch  the  properties  of  Householder  transformations  here. 

For  computational  details  see  [3,6,7 J . 

A Householder  transformation  is  a symmetric,  orthogonal  matrix  of 

the  form 


For  any  nonzero  vector  a there  is  a Householder  transformation  H such 

a 


that 

(3.1)  H a = T il  ail 2 er 


100 


- 6 - 


where  e^  = (1,0,. .. ,0)^.  In  fact  we  may  take  u = e^  ± a,  where,  for  reasons 
of  numerical  stability,  the  sign  is  chosen  to  make  ± aj  >0  . It  follows  from 

(3.1)  that  the  first  column  of  H is  the  vector  a scaled  to  have  norm  one, 

d 

while  the  remaining  columns  form  an  orthonormal  basis  for  the  space  orthogonal 
to  A. 

Householder  transformations  may  be  stored  and  manipulated  easily. 

For  their  storage,  they  only  require  space  to  retain  the  vector  u.  Moreover, 
given  u,  the  product  Ha  can  be  calculated  for  any  vector  a in  the  form 

Ha  = a - ^-4-  u, 

Hull  2 

without  the  necessity  of  forming  H explicitly. 

The  decomposition  we  will  estimate  from  is  described  in  the 
following  theorem. 

Theorem  3.1.  Let  A be  a matrix  of  order  n and  let  J be  a 
permutation  matrix.  Then  there  is  an  orthogonal  matrix  Q and  an  upper  tri- 
angular matrix  R such  that 

A = QRJ. 

Moreover  J can  be  chosen  so  that 

(3.2)  r^  > ^ r.jj  (k=l , . . . ,n  j j=k, . . . ,n) . 

Proof.  We  will  show  how  J can  be  chosen  so  that  (3.2)  is  satisfied, 
the  proof  for  an  arbitrary  fixed  J being  similar.  Let  the  columns  of  A 
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be  rearranged  so  that  one  of  largest  2-norm  is  first.  This  amounts  to 
choosing  a permutation  matrix  0^  so  that 

Aj  = Aj{  = (aj1),a^1),...,aj1)) 


satisfies 


(3.3) 


41  "2  - Haj  H2  (j  = 2’3’ 


Let  be  a Householder  transformation  such  that 

H a ( 1 ) e 

Vl  rller 


Then  HjAj  has  the  form 


Vi- 


rn  rI 

0 A0 


where  r1  is  an  (n-l)-vector  and  A is  of  order  n-1.  By  an  induction 
hypothesis,  we  may  assume  that  there  is  an  orthogonal  matrix  Q2  and  a 
permutation  matrix  J2  such  that 

^2^232  = ^2’ 


where  the  elements  of  R2  satisfy  (3.2).  If  we  set 


■«.c 


1 0 


iO  J, 


, • ? 
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then 


qtaot  * 


rll 
I o 


R2  ' 


= R 


The  matrix  R is  clearly  upper  triangular.  To  show  that  it  satis- 
fies (3.2),  note  first  that  since  R2  also  satisfies  (3.2),  it  is  only  nec- 
essary to  treat  the  case  k = 1.  Now  r^  = ila^ll.  Moreover  if  the  j-th 
element  of  rjj2  corresponds  to  the  z-th  column  of  A^,  then 
||a|^||  2 = rij  + r2j  +--'+  rjj-  The  result  now  follows  from  (3.3).  a 

When  J is  fixed,  the  above  proof  is  essentially  a specification  for 
an  algorithm  to  factor  AJT  into  the  product  QR  of  an  orthogonal  matrix  and 
an  upper  triangular  matrix.  If  J=I,  this  is  called  the  QR  factorization 
of  A.  For  A nonsingular,  this  factorization  is  unique  if  the  diagonal  ele- 
ments of  R are  required  to  be  positive. 

The  process,  sketched  in  the  proof,  by  which  columns  of  maximal  norms 
are  moved  to  the  front  of  the  current  matrix  is  called  column  pivoting.  We 
will  call  the  resulting  decomposition  the  QRJ  factorization  of  A.  Unless 
there  are  ties  among  norms  of  columns  during  the  reduction,  it  is  unique  up  to 
the  signs  of  the  diagonal  elements  of  R.  Because  of  the  condition  (3.2), 
the  matrix  R tends  to  show  a graded  structure  with  the  elements  becoming 
smaller  as  one  procedes  toward  the  lower  right  hand  corner.  In  particular, 
if  rank(A)=r  then  r. ,=0  for  j > i > r.  Near  singularity  of  A will 

’ J 

manifest  itself  in  a large  ratio  of  r^  to  rnr).  This  suggests  that  the 
number  r^/rnn  may  provide  an  estimate  of  <2(A);  i.e., 


It  is  this  estimate  of  < ^ that  we  shall  study  in  this  paper. 

4.  Random  Orthogonal  Matrices. 

The  matrices  used  in  our  study  are  generated  by  transforming  diagonal 
matrices  by  random  orthogonal  matrices.  Accordingly,  in  this  section  we  de- 
scribe an  economical  method  for  generating  random  orthogonal  matrices. 

Just  as  there  is  a natural  uniform  distribution  for  real  numbers,  so 
is  there  a natural  distribution  associated  with  the  space  of  orthogonal  matrices. 
Let  denote  the  set  of  orthogonal  matrices  of  order  n.  With  the  spectral 
norm  as  a metric,  & is  a compact  topological  group.  Hence  there  is  a unique 
normalized,  left-translation  invariant  measure  p for  $n;  that  is,  p 
satisfies  p(f>n)  = 1 and 

u(H*)  = p(*) 

for  any  measurable  % c.  & and  any  He  6^15}.  It  can  be  shown  that  p is 
also  right  translation  invariant.  In  the  sequel  all  functions  and  distribu- 
tions will  be  assumed  to  be  measurable  with  respect  to  p and  all  integrations 
will  be  performed  with  respect  to  p. 

Our  method  for  generating  random  orthogonal  matrices  is  based  on  the 
following  observation.  Let  X be  an  n x n matrix  whose  elements  are  inde- 
pendently normally  distributed  with  mean  zero  and  common  variance  (we 
shall  abbreviate  this  statement  by  saying  that  X is  distributed  NI(0,cj2)). 


mu  JJi  - . _■■  ■ 
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Let  X = QR  be  the  QR  factorization  of  X,  normalized  so  that  the  diagonal 
elements  of  R are  positive.  Then  Q is  distributed  over  fr  according  to 
y . 

This  fact  follows  from  the  invariance  of  the  distribution  of  X 
under  orthogonal  transformations.  Specifically,  let  H be  a fixed  orthogonal 
transformation.  Then  HX  has  the  same  distribution  as  X.  Since  HQ  is 
orthogonal  and  HX  = (HQ)R,  the  random  matrix  HQ  is  the  orthogonal  part  of 
the  QR  factorization  of  HX.  It  follows  that  HQ  has  the  same  distribu- 
tion as  Q;  i.e.  the  distribution  of  Q on  5 is  invariant  under  left 
translations.  Hence  Q has  the  distribution  y. 

All  this  suggest  that  random  orthogonal  matrices  be  generated  by 
generating  a matrix  X from  NI(0,o?)  and  computing  its  QR  factorization. 
However  this  ocess  requires  0(n3)  operations.  We  shall  show  that  we  can 
instead  express  a random  orthogonal  matrix  as  a product  of  Householder  trans- 
formations that  can  be  generated  in  0(n2)  operations.  We  begin  our  devel- 

2 

opment  with  an  elementary  lemma,  in  which  we  extend  our  notation  NI(0,a  ) 
to  refer  to  vectors  as  well  as  matrices. 

2 

Lemma  4.1.  Let  the  random  n-vector  x be  distributed  NI(0,a  ) and 

let  H be  a random  matrix  on  that  is  independent  of  x.  Then  y = Hx 
2 

is  NI(0,o  ) and  independent  of  H. 

Proof.  Let  f and  g denote  the  distributions  of  x and  H 
respectively,  and  let  h(y,H)  denote  the  joint  distribution  of  y and  H. 
Let  and  H be  measurable  sets.  Then 


11  - 


(4.1) 

Now  for  fixed  H 


f fh(y,H)  dydH  = Pr  [ (y ,H)  fcijx*]. 

* 3f 


Pr[y^/lH]  = Pr  [x  e HTi/  I H] 

= Pr  [ xt  yi  = ( f(x)dx. 

H 

Then  next  to  last  equality  follows  from  the  orthogonal  invariance  of  the 
distrib-tion  of  x and  the  independence  of  x from  H.  It  follows  from 
(4.1)  and  (4.2)  that 


( h(y,H)  dydH  = ( g(H)  Pr  l y fc  ^ \ h] 

n -HI 


dH 


f(x)  dxdH. 


Hence  h(y,H)  = f(y)g(H)  almost  everywhere,  which  implies  the  result.  ■ 

The  computational  algorithm  is  based  on  the  following  theorem.  Its 

complete  statement  is  due  to  Garrett  Birkhof,  although  the  author  discovered 

parts  independently  in  the  course  of  deriving  the  algorithm  to  follow. 

Theorem  4.2.  Let  X be  distributed  NI(0,c^)  and  let  X = QR  be 

the  QR  factorization  of  X normalized  so  that  the  diagonal  elements  of  R 

are  positive.  Then  Q has  the  distribution  p on  & . The  elements  of  R 

are  independent  of  Q and  of  each  other.  In  particular  the  super  diagonal 

2 2 2 2 

elements  of  R are  NI(0,a  ),  while  r^./o  has  x distribution  with 
n-i+1  degrees  of  freedom. 
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Proof.  We  consider  a variant  of  the  first  step  of  the  reduction  to 
triangular  form  sketched  in  Section  3.  From  the  first  column  x1  of  X,  we 
determine  an  orthogonal  matrix  (which  is  ± a Householder  transformation) 
such  that  HjXj  = il x ^11  e ^ . Then  HjX  has  the  form 

/rU  r12  rln  l 


(4.3) 


HjX  = (HjXj, 


H1x2,...H1xn 


\ _ 
> 


,(2) 


.0) 


0.9  2 

Now  Hj  depends  only  on  Xj/liXjH,  while  r^  = il  x^li  . Since  Xj,  is  Nl(0,a  I), 

2 2 

it  is  easily  seen  that  and  r^  are  independent  and  moreover  r^/a 
2 

has  x distribution  with  n degrees  of  freedom. 

By  Lemma  4.1  the  vectors 


ru  i 

= H.x,  (j  = 2 n) 

xf>)  J 

o 

are  independent  of  each  other  and  H1  and  are  NI(0,o  ).  It  follows  that 
(r^.- • • »^jn)  and  X2  = (x2  ,...,x^  ')  are  independent  and  NI(0,o  ). 

If  the  reduction  is  continued  with  X2>  there  results  a sequence  of 
orthogonal  matrices,  such  that 

Hn  HjX-R. 

where  R is  upper  triangular.  By  the  arguments  of  the  last  paragraph,  the 
elements  of  R are  independent  of  H^,. . . , Hn  and  have  the  distributions 
claimed  in  the  theorem  statement.  The  matrix  Q = is  the  orthogonal 

part  of  the  QR  factorization  of  X and  by  the  discussion  above  must  have 
the  distribution  p on  V . ■ 


I 
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The  computation  of  Hj,  H2,...,  HR  as  described  in  the  proof  of 
the  theorem  involves  the  reduction  of  the  entire  matrix  X and  consequently 
entails  O(n^)  operations.  However,  note  that  depends  only  on  x^, 

which  can  be  generated  without  generating  the  other  columns  of  X.  Likewise, 

(2)  2 
depends  only  on  x2  , which  is  simply  a vector  distributed  NI(0,o  ). 

lo) 

Thus  it  is  unnecessary  to  generate  x2  and  compute  HjX2-  Rather,  x2 
can  be  generated  directly,  and  H2  computed  from  it  — and  so  on  for  the 
rest  of  the  H's. 

There  is  a minor  annoyance  to  remove.  The  H's  in  the  proof  of 
Theorem  4.2  are  not  Householder  transformations;  rather  they  are  Householder 
transformations  scaled  by  ±1  to  make  the  diagonal  elements  of  R positive. 
However,  it  is  easily  seen  that  the  factor  is  independent  of  the  Householder 
transformation,  since  it  depends  only  on  the  sign  of  the  first  element  in  the 
generating  vector  xj^.  Hence  the  matrix  of  the  theorem  differs  from  the 
product  of  the  Householder  transformations  by  a factor  of  ±1  that  is  indepen- 
dent of  the  product.  Such  a factor  does  not  change  the  distribution;  hence 
in  generating  Q we  need  work  only  with  the  Householder  transformations. 

We  summarize  the  above  considerations  in  the  following  theorem. 
Theorem  4.3.  Let  the  independent  vectors  Xj.Xg xn  be  dis- 
tributed NI(0,a2)  over  !Rn,  Rn-1,. ...  R1.  For  j = 1,2 n,  let  fl. 

XJ 

be  the  Householder  transformation  that  reduces  x.  to  *Hx.lle.  cf.  (3.1). 

J J • 

Let  Hj  = diag(Ij_ j.flj) . Then  the  product  Q = H^  H2...  Hn  is  a random 
orthogonal  matrix  distributed  according  to  the  Haar  measure  over  ^n- 


* 
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The  procedure  implied  by  Theorem  4.1  is  quite  economical.  The 
Householder  transformations  require  0 ( n ) operations  to  generate,  and  it 
is  probable  that  most  of  the  work  will  be  spent  generating  pseudo-random 
normal  deviates  for  the  vectors  x..  If  Q is  required  explicitly,  the 

*J 

transformations  must  be  multiplied  together,  which  is  an  0(n  ) process. 

However,  in  many  applications  this  will  be  unnecessary.  For  example,  if 

the  product  Qx  is  needed,  it  can  be  formed  by  successively  multiplying 

H . H,  into  x,  an  0(n  ) process, 

n n- 1 i. 

5 . The  Experiment  and  Results . 

In  order  to  test  the  condition  estimaters  described  in  Sections  2 
and  3,  we  generated  random  test  matrices  as  follows.  Given  a diagonal  matrix 


£ — di ag  0^) 

satisfying  > a ^ >...z  > o,  random  orthogonal  matrices  U and  V 

were  generated  and  A was  computed  as 

A = UeVT . 


Note  that  < ^ (A)  = so  that  the  condition  number  of  A can  be  controlled 

by  adjusting  and  Further  replications  of  A having  the  same  condi- 

tion number  can  be  obtained  by  holding  I fixed  and  generating  new  U and  V. 
Three  factors  were  varied  in  the  experiment. 

1.  The  order  of  the  matrix.  We  took  n = 5,10,25,50  . 


2.  The  condition  number.  We  took 


10,102,104,106. 


The  distribution  of  the  singular  values  a.. . We  considered 


two  distributions.  In  the  first  1 


n-1' 


-1 
c2  • 


i — — *— -*p 


3. 


This  represents  the  case  where  the  ill -conditioning  of  A 
is  due  to  a sharp  break  in  the  singular  values.  In  the  second 
distribution  we  took  = 1,  an  = k^1,  and  a2^°l  = a3^a2  = • • • 

= 0n/°n_i*  In  this  case  the  singular  values  decay  exponentially 
from  Oj  to  Op,  and  it  is  impossible  to  point  to  a single  one 
that  causes  the  ill-conditioning. 

For  each  case,  corresponding  to  a choice  of  the  three  factors, 
twenty-five  matrices  were  generated  by  varying  U and  V randomly  as  de- 
scribed above.  Estimates  Kj  and  were  calculated  as  described  in 
Sections  2 and  3 via  UNPACK  codes  SGECO  and  SQRDC,  and  the  ratios 

and  <2^k2  recorded*  The  number  was  given  as  one  of  the  factors  of 
the  case.  The  number  was  computed  by  inverting  A.  Care  was  taken 
throughout  the  experiment  to  insure  that  rounding  error  played  an  insignif- 
icant part. 

The  raw  data  is  too  voluminous  to  present  here.  It  is  summarized  in 
the  histograms  in  the  appendix  to  this  paper  (the  scale  is  xlO).  Table  1 
gives  medians  of  the  ratios  estimated  from  the  histograms. 

The  table  shows  that  for  the  classes  of  matrices  treated  here,  both 
techniques  do  a reasonable  job  of  estimating  the  condition  number.  The 
algorithm  in  SGECO  generally  outperforms  the  QRJ  factorization,  as  might  be 
expected  since  SGECO  was  specifically  designed  as  a condition  estimater.  In 
all  cases  the  quality  of  the  estimate  deterioriates  with  increasing  n,  but 
not  drastically  so.  In  the  entire  experiment  no  ratio  less  than  0.1  was 


observed. 


Table  1 

Median  of  the  Ratios 


exponential  decay  exponential  decay 
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The  distribution  of  the  singular  values  has  a marked  effect  on  the 
condition  estimaters.  Both  do  better  when  there  is  a sharp  break.  For  this 
distribution,  the  condition  number  itself  has  little  effect  on  the  quality 
of  the  estimate.  For  the  exponential  decay,  the  two  estimaters  behave  dif- 
ferently. SGECO  improves  as  the  condition  number  increases,  while  the  QRJ 
factorization  deteriorates. 

It  is  risky  to  make  sweeping  generalizations  on  the  basis  of  experi- 
ments as  limited  as  these.  However,  to  the  extent  that  they  are  applicable 
they  suggest  that  the  techniques  will  be  useful  when  an  order  of  magnitude 
estimate  of  the  condition  number  is  desired.  This  is  reinforced  by  the  fact 
that  to  date  no  one  has  succeeded  in  making  the  algorithms  fail  dramatically. 
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APPENDIX 


Histograms  of  Condition 
Number  Ratios 
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